Theory of doped Mott insulators: duality between pairing and magnetism 
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By bosonizing the electronic t-J model exactly on any two-dimensional (2D) lattices, and inte- 
grating out the gauge fluctuations combined to slave particles beyond mean fields, we get a theory in 
terms of physical Cooper pair and spin condensates. In the sense of mutual Berry phase they turns 
out to be dual to each other. The mutual-duality is the missing key in the resonant-valance-bond 
idea[l] to work as a paradigm of doped 2D Mott insulators. We argue that essential aspects of 
high-Tc phenomenology find natural solutions in the theory. We also provide interesting predictions 
for systems on hexagonal lattices. 



/. Introduction: Ever since Anderson's proposal that 
superconductivity in high-Tc copper oxides may be a real- 
ization of resonant- valence-bond (RVB) physics, [1] solv- 
ing the t-J model, [2, 3] a standard electron model for 
systems with strong on-site Coulomb repulsion, has been 
not only a daydream but also a nightmare in the commu- 
nity for more than a decade. The good news is that in- 
deed a lot of high-Tc phenomenology may be interpreted 
in one way or another using the t-J model. The bad news 
is a consensus to all phenomenology is not reached even 
just in the underdoped region. It is therefore academ- 
ically desirable to have a theoretical paradigm (in the 
infrared limit), just as Luttinger-liquid theory in one- 
dimensional (ID) systems, Landau-Fermi liquid theory 
for metals and quantum Hall liquid for semiconductor 
wells under high magnetic fields. On the theory side 
such a paradigm helps to answer the question: Does t-J 
model describe high-T c anyway? If it does, what makes 
high-Tc unique? Inspired by the fact that Luttinger liq- 
uid theory can be easily brought out by bosonization in 
ID and recent successes of the phase-string theory, [3] in 
this Letter we tackle the challenge raised by Anderson 
to generalize the Luttinger liquid to 2D strongly corre- 
lated electron systems. [4] We bosonize the electronic t-J 
model exactly [see Eqs.(9) and (14)]. The mean field 
theory integrates spin ordering and pairing thanks to 
the "all-boson" hamiltonian. We further integrate out 
gauge fluctuations beyond mean field and obtain a the- 
ory in terms of physical Cooper pair and spin magnetic 
condensates alone [see Eqs.(22), (23), and (25)]. They 
turn out to be dual to each other in the sense that 
Cooper pair condensate (spin condensate) sees the classi- 
cal part of magnetic moments (Cooper pairs) as 27r-flux 
bundles. [5]. The mutual-duality between superconduc- 
tivity and magnetism is clearly what makes doped 2D 
Mott insulators unique, and is just the missing key in 
the RVB theory. Because of the intervening magnetism, 
it seems unnecessary to prepare a quantum spin liquid 
in parent Mott insulators to adiabatically connect su- 
perconductivity. Other far-reaching consequences of the 
theory are detailed in the closing section after warming 
up with the notations and technical expositions, which 



could be skipped for anxious readers. 

II. Bosonization of the t-J model: The electronic 
t-J model in the so-called slave-boson representation can 
be written as 
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In the above hi (fia) is hardcore holon (fermionic o- 
spinon annihilation operators). We have written Hj in 
the hopping and singlet pairing channels. They help iden- 
tify the saddle point of the theory in any dimension and 
for any lattice geometry. The physical Hilbert space is a 
subspace of the Fock space subject to 



hlh + J2flfi<r = -i- 
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The slave-boson hamiltonian is invariant under the local 
U(l) gauge transform 



ficrexpiiOi), hi — > h t exp(i9i). 



(6) 



This symmetry guarantees the conservation of the local 
gauge charge h\hi+Y^, a f tafia and the associated current. 

While not necessarily unique, we define the following 
Jordan- Wigner transformation 
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were N a is the Klein factor, bi„ hardcore boson anni- 
hilation operator, 9u is a multi-valued function defined 
by 
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n ia — bj a bicr and finally n a is the average occupation of 
CT-spinons. For systems with anti-ferromagnetic exchange 
it is reasonable to set n a = (1 — x)/2 with x the holon 
concentration. The purpose of including n a is to make 
the gauge field attached to the 6 CT -bosons vanish at the 
mean field level. 

By simple substitution and taking into account of the 
hardcore nature of the bosons as well as the occupation 
constraint, Hj can be rewritten as 
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Here = h\hi. The bond phase J7»,i±5 = =F(l~ a; ) 7r /2 for 
5 = x and y on square lattices, and for 5 = x and y/3y/2i 
x/2 on hexagonal lattices. From the above we see that 
holons appear as n-flux bundles as viewed by spinnons. 
This effect is cancelled by the background field due to x 
at the mean field level, but should be kept beyond the 
mean field level due to its topological nature. Certainly 
the above definition of /Cy is unique up to arbitrary spin- 
independent link phase factors, 
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for any oy, a reflection of the local U(l) gauge symmetry. 
We choose not to present the pseudo-pairing operator at 
this stage. As will be clear, pairing is an inevitable mode 
in the theory. 

Now we can proceed to transform the kinetic part H t . 
By substitution and manipulating in terms of , we get 
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where we have used the following identity 
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Therefore z-direction magnetic moments appear as Tr-flux 
bundles as viewed by holons. Again the above definition 
of K. h is unique up to the same U(l) gauge transform 
applied to JQj . 

Summarizing, H = H t + Hj — fi^^l — n\) with the 
new forms of H t and Hj in Eqs.(9) and (14) discussed 



above forms another rigorous fully-bosonized represen- 
tation of the t-J model. The holons and spinnons are 
mutually frustrated as reflected by the gauge fields a h 
and a m . We note that a similar notion was pointed out 
in the phase string theory designed particularly for bipar- 
tite lattices. [3] The present theory applies for any lattice 
geometries. 

III. Mean field theory: The interacting part 
of the hamiltonian can be rewritten as H t + Hj = 
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with Hij = Kij 



y/C^ . This suggests the MF effective hamiltonians for 
holons and spinons, respectively, 

H h MF = -i]T(x^+h.c.)+ M ]>>f, (18) 
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where we suppressed the residual interaction for brevity. 
Here \ij IS the mean field hopping amplitude \ij = (Hij)- 
The Lagrangian multipliers [i and A enforces (h\hi) = x 
and ^2i a (b\ a bia) = 1 — x, respectively. Let us assume 
zero average magnetization in the z-direction so that both 
holons and spinons see zero flux in the MF theory. An 
obvious mean field solution is \ij — s \x\ with s — sign(t) 
so that all three species of bosons can condense to their 
own non- degenerate ground states: hi — y/x, and bi a = 

yj ^~2^ cxp[— iasQ ■ ri\ with Q = Qo(l — sx). Here Q = 

(it, 7r) and 27r(l/3, 1/y/S) on a square and hexagonal lat- 
tice, respectively. By self consistency, |x| = l—x+2\t\x/J 
and (1 — x) sin[(l — sx/2)n/3] + 2\t\x/J on square and 
hexagonal lattices, respectively. Q would be the mag- 
netic ordering wave vector as seen from (S^~) = (b^hi) = 
^(1 — x) exp(isQ • Ti). The s-dependence in Q on square 
lattices is trivial because Qo(l ± x) are identical, but is 
nontrivial on hexagonal lattices. One would also thought 
of spin pairing b^bji in the condensates. The bosonized 
MF state therefore integrates hopping, magnetism and 
spin pairing. However, we must remind ourselves that 
the mean field breaks the local U(l) gauge symmetry. 
It is a must to integrate out gauge fluctuations beyond 
mean field theory to talk about physical entities. This is 
the subject of the gauge theory in the next section. 

IV. Integration over gauge fluctuations: The Eu- 
clidian action S = J Ld 2 xdr that describes the long- 
wave-length low-energy fluctuations beyond the above 
MF state is as follows. The Lagrangian density L = 
L h+J2a L a +L cs +ia with L h , L a and L cs describe, re- 
spectively, holons, (7-spinons and the Chern-Simons-likc 
flux-attaching reflecting the mutual frustration between 
holons and spinons. The holon part is given by 
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In the above ph = x is the average holon concentration, 
Sph is the holon density fluctuation, j/j is the holon cur- 
rent density, K h oc 2tx\ is the bare holon phase stiffness 
(this is best seen by integrating out j/j), ^ is defined by 
h = \h\<f>h, (ao,a) is the auxiliary U(l) gauge field, a[f 
is a Lagrangian multiplier for flux attaching to frustrate 
spinons, a m is the gauge field due to the frustration from 
spinons, (A ,A) is the physical electro-magnetic (EM) 
gauge field coupled to the positively-charged holons, and 
finally Uh is an effective short-range repulsion between 
holons. Similarly, the tr-spinon part L a is given by 

L a = i(p a +Sp a )((t)* a — (j) a - a Q ) -iaSp a a r Q +^-dpl 

+ i U ■ (C-0. - a - ,a ft ) + JL,£. (21) 

In the above a™ is a Lagrangian multiplier that at- 
tach flux to the z-direction magnetization that frustrates 
holons, Kb oc Jx(l — x)/2 is the bare spinon phase stiff- 



ness, and is defined by b„ = \b a \ exp(— iaQ ■ r/2)(f> a 
which describes fluctuations of the cr-boson beyond the 
condensate momentum. The last piece Lcs is given by 

Lcs = ~4(y x a"), + -a™(V x a m ) z . (22) 
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In the above the subscript z means the z-component of 
the object. Let us observe that integration over a , Oq 
and a™ using the total action L indeed enforces occupa- 
tion constraint and attaches 7r-flux to holon and magne- 
tization fluctuations, respectively. 

Given the above form of the action, we can integrate 
out the auxiliary U(l) gauge field (do, a) exactly ,[6] 
yielding Sp h + o><x = and ih + E<x 3° = °- Tnis 
enables us to get rid of slave particles in favor of a U(l) 
gauge-neutral theory L = L c + L rn + Lcs as follows. 

First, L c describes the charge condensate (CC) of the 
theory, 
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Here we have used the identity (f> a d t (f>a ~ 2<p* h d t 4>h = 4>* c dt^ c - Clearly 4> c is the phase factor of bib±h*h*, a spin- 
charge recombined Cooper pair. Moreover, p c and j c are just the Cooper pair density and spatial current, respectively, 
K c is the effective zero-temperature superfluid density, and finally u c is the effective local charge repulsion. At low 
doping K c oc x. Interestingly, Aq in the first two terms of L c simply reminds us that the effective charge density is 
1 — 2p c — x, but the Cooper pairs response to EM field as 2e-carriers with pair density p c = (1 — x)/2. This proves 
that pairing does not imply bounding of holes at all (they are not the normal modes in the theory), a conclusion 
supported by variational Monte Carlo studies of the t-J+Coulomb model. [7] Moreover, the Cooper pairs also view a 
localized z-direction magnetic moment as a 2n-flux bundle as described by 2& m . 
Second, L m describes the spin magnetic condensate (MC) of the theory, 

8 V u 1 

L m = iSp^Y^m - 2<C) + j m • (CjAn - 2a ft ) + -fSp 2 m + —f m , (25) 

S Pm = - ^aSpv, (f) m = <P]<P\i jm = ^ Crjcr, K m = -K b , u rn = 2u b . (26) 



Clearly, 8p m , <j) m and j m are the z-direction magne- 
tization S z , the phase factor of S~ and the spin cur- 
rent (beyond the magnetic ordering wave vector Q). The 
gauge field 2a. h describes the effect of Cooper pair density 
fluctuations on the spin dynamics. A localized Cooper 
pair would appear as a 27r-flux bundle to S~ . On the 
other hand, MC condensate is not coupled to physical 
EM gauge fields and is therefore charge-neutral, as would 
be desirable. 

Except for the first two terms in L c (which gives the 
average density of Cooper pairs and the charge carrier 



density) the two condensates CC and MC are mutually 
dual to each other. This is the crucial part of the theory 
that makes a difference to the ID counterpart in which 
spin-charge separation is complete. 

V. Discussions: The following discussion at zero 
temperature relies on the notion of mutual-duality iden- 
tified above and the self-duality for a separate conden- 
sate alone. [5, 8] We specialize to the context of high-T c 
phenomenology, [9] but the basic physics applies to all 
doped Mott insulators. 1) At zero doping we have a quan- 
tum antifcrromagnetism described by the MC alone. The 



4 



gauge field a h = since there is no charge fluctuations 
at all. 2) At low doping CC is disordered by quantum 
(or thermal) fluctuations due to the small phase stiffness. 
The charges are localized and act as vortices to MC by 
mutual-duality. Depending on whether MC is of type- 
I or type-II "superconductor" , these charges may either 
phase separate or form a Wigner crystal, respectively. 
They correspond to the intermediate state and mixed 
state in usual superconductors, and would suggest two 
types of Mott insulators as recently advocated. [10] This 
is a useful concept to understand the absence or pres- 
ence of electronic inhomogeneities in different families of 
high-Tc copper oxides. [11] Localized charges can make 
the spin order glassy. Another interesting possibility is 
the charge stripe order in cuprates. Since the localized 
charges appear as vortices to MC, the fact that stripe is 
an anti-phase boundary to the underlying magnetic or- 
der [12] seems very natural. 3) With increasing doping, 
the charges are more and more mobile and eventually 
CC shows up above some lower critical doping level. 3a) 
As we envisioned above, although the doped holes do 
not pair at all, the Cooper pairs response to EM fields 
as charge-2e carriers. This determines the correct flux 
quantum as hc/2e. 3b) According to the composition 
rule of supcrfluid density K c in Eqs.(24), for small x it is 
dominated by the smaller component Kh , and is therefore 
roughly linear in x. 4) Consider ordered CC and disor- 
dered MC, a superconductor with no spin order. The 
vortices of MC must condense by self-duality. [5] Let us 
call this new condensate as M'C. Just as CC does, M'C 
also views localized magnetization S z as 27r-flux bun- 
dles, i.e., the next generation vortex excitations. There- 
fore CC and M'C can happily coexist. This opens up 
a handful of interesting possibilities. 4a) Vortices un- 
der an applied magnetic field in CC may capture local 
spin magnetic moments, forming non-topological vortices 
with zero winding number in CC to lower the superflow 
kinetic energy. Since winding is a necessary condition to 
the zero energy electron state in a vortex core of d-wave 
superconductor, [13] non-topological vortices may resolve 
the puzzle that no vortex-core zero-bias tunnelling con- 
ductance peak has ever been observed in cuprates. [14] 
4b) A vortex in CC may also capture many spin mo- 
ments (favorably ordered with staggered polarization), 
and this would lead to the checkerboard pattern in the 
local density of states near a vortex core in cuprates. [15] 
4c) Moreover, since the captured spins are vortices to 
M'C by self-duality, it is favorable to recover MC. This 
may explain the enhancement of anti-ferromagnetic cor- 
relation in the vortex state of cuprates. [16] 4d) The spin 
excitations (Sp m ,j m ) appear as the "EM" fields to M'C 
by self-duality. [5] Such "photons" are gapped up to the 
plasma frequency of CC+M'C state that scales with the 
physical superfluid density. In view of our starting mag- 
netic vector Q, this gap actually corresponds to the in- 
elastic neutron resonance energy found in cuprates. [17] 



The doping dependent incommensurability in Q is obvi- 
ous in the theory. 4e) On the other hand, assuming the 
absence of long-range interactions in either CC and/or 
M'C, longitudinal phase fluctuations in the condensates 
are phonon-likc, and via rc-fermionization of the emerg- 
ing low energy field theory[18] they are nothing but the 
Dirac-like nodal quasi-particles in a d-wave superconduc- 
tor. 4f) Topologically de-confined but spatially localized 
spin moments, such as caused by Zn, induces vortices in 
CC at zero magnetic field by mutual-duality. Therefore 
the critical Zn concentration to kill CC would be given 
by H c2 /^o- (The spins do not have to point in the same 
direction.) This seems to agree with the so-called "Swiss 
cheese" model of Zn in cuprates. [19] More detailed dis- 
cussions are in progress. 
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